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Abstract 



A search for physics beyond the standard model involving events with one or more 
photons, jets, and missing transverse energy has been performed by the CMS ex- 
periment. The data sample corresponds to an integrated luminosity of 4.93 fb _1 of 
proton-proton collisions at y/s = 7 TeV, produced at the Large Hadron Collider. No 
excess of events with large missing transverse energy is observed beyond expecta- 
tions from standard model processes, and upper limits on the signal production cross 
sections for new physics processes are set at the 95% confidence level. The results of 
this search are interpreted in the context of three models of new physics: a general 
model of gauge-mediated supersymmetry breaking, Simplified Models, and a theory 
involving universal extra dimensions. In the absence of evidence for new physics, 
exclusion regions are derived in the parameter spaces of the respective models. 
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1 Introduction 

The standard model (SM) of particle physics is a very successful theory describing existing ex- 
perimental data. However it is not expected to describe physics up to the Planck scale, because 
of the extreme fine tuning required to control particle masses (hierarchy problem) HHHl, nor 
does it provide an explanation for dark matter. These issues with the SM motivate a broad pro- 
gram of searches for physics beyond the SM. Among the theories proposing physics beyond 
the SM, supersymmetry (SUSY) is of particular interest as it resolves these problems by intro- 
ducing a symmetry between fermions and bosons resulting in a superpartner (sparticle) for 
each SM particle with identical quantum numbers except spin. Since no sparticles have been 
found yet, SUSY must be a broken symmetry with the masses of the supersymmetric particles 
being heavier than their SM partners. The version of supersymmetry based on general gauge- 
mediated (GGM) SUSY breaking E HlOU is of particular theoretical interest for new physics as it 
not only stabilizes the mass of the SM Higgs boson and drives the grand unification of forces, 
but also avoids the large flavor-changing neutral currents that trouble other SUSY-breaking 
scenarios. Another extension to the SM is the theory of universal extra dimensions (UED) (3|, 
which predicts additional compactified dimensions beyond the regular four space-time dimen- 
sions of the SM. These extra dimensions (ED), which are accessible to standard model fields, 
could allow gauge coupling unification and provide new mechanisms for the generation of 
fermion mass hierarchies. 

This paper describes a search for events with two signatures containing photons, which may 
indicate new-physics processes in a variety of theoretical scenarios including GGM supersym- 
metry and UED. Final states with photons are experimentally interesting as photons can be 
identified with relatively high purity and efficiency with the Compact Muon Solenoid (CMS) 
detector. The first signature studied consists of at least one isolated photon with high en- 
ergy measured in the plane transverse to the beam direction (Ej), at least two hadronic jets, 
and large missing transverse energy (E™ 1SS ). The second signature is characterized by at least 
two isolated photons with high Ej, at least one jet, and large E™ ss . This search is based on a 
data sample recorded with the CMS experiment corresponding to an integrated luminosity of 
4.93 ± 0.11 fbr 1 of pp collisions at = 7TeV produced at the Large Hadron Collider (LHC). 

The organization of this paper is as follows. This introductory section is followed in Section [2] 
by a discussion of the theoretical framework used for the interpretation of this search, and then 
in Section [3] by a description of the CMS detector. The event selection criteria are detailed in 
Section |4] and the description of the simulated samples is given in Section [5] The methodology 
to estimate backgrounds is explained in Section|6j Sections |7]and|8] discuss details of the single- 
photon and diphoton analyses including the experimental results. Section |9.1| expresses the 
search results in terms of exclusion regions in the context of the GGM SUSY scenario, while 



in Section 9.2 and Section 9.3 the results are interpreted in the context of a final state driven 
"simplified" model, and universal extra dimensions, respectively. Conclusions are stated in 
Section M 



2 Theoretical Framework 



The result of this search is interpreted in the context of three models of new physics. We discuss 
in this section the theoretical framework on which these models are based. 
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2 Theoretical Framework 




2.1 General Gauge-Mediated Supersymmetry Breaking 

The first model is a gauge-mediated SUSY scenario ||TT1 - [T3| in which the gravitino (G) is the 
lightest SUSY particle (LSP) and the lightest neutralino (x^) is the next-to-lightest SUSY particle 
(NLSP). The gravitino escapes detection, leading to E™ 1SS in the event. The neutralino in the 
GGM models that we consider consists predominantly of either the bino, the superpartner of 
the U(l) gauge field, or the wino, the superpartner of the SU(2) gauge fields. Assuming that 
R parity |T4| is conserved, strongly-interacting SUSY particles are pair-produced at the LHC. 
Their decay chain includes one or more quarks and gluons, as well as the neutralino NLSP, 
which in turn decays into a gravitino and a photon or a Z boson. Figure [T] shows several 
diagrams of possible GGM processes that result in a single-photon or diphoton final state, in 
squark and gluino pair production processes. If the NLSP is bino-like, its branching fraction to 
a photon and gravitino is expected to be large, leading to an enhancement of the diphoton final 
state (see Fig.[I]top). If the NLSP is wino-like, its branching fraction to a photon and gravitino is 
reduced, leading to a relative enhancement of the single-photon final state (see Fig. [l]bottom). 
Therefore we perform searches in both the single-photon and diphoton final states in order to 
be sensitive to models with different NLSP composition. 

Table [T] provides examples of such GGM decay chains leading to photons in the final state. 
The table is divided horizontally between single-photon and diphoton final states. The vertical 
direction differentiates between bino-NLSP and wino-co-NLSP cases. The number of jets pro- 
duced in the cascades can vary depending on whether gluinos or squarks are produced, and 
the species of quarks in the final state. This search is also sensitive to the scenario in which the 
NLSP is a pure wino. In that case, the lightest chargino (x~\) is also a wino, and the chargino- 
neutralino mass difference is too small for one to decay into the other, resulting in the chargino 
to decay directly into a gravitino and a W boson (see Fig.JT]). In this analysis we do not veto on 
the presence of isolated leptons since in the wino co-NLSP case we seek to detect the neutralino 
decays into Z bosons and chargino decays into W ± , both of which decay chains can result in 
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Table 1: Examples of GGM cascades leading to the topologies of a single photon or diphotons 
in the final state. 



NLSP type 


7 + 2 jets + E™ 1SS 


77 + jet + EJjF ss 


Bino-like 


jets + x ^ -> jets + 7 + Z + GG 


jets + j ets + 77 + GG 


Wino-like 


jets + x ^ -> jets + 7 + Z + GG 
jets + x^Xi ~» j' ets + 7 + W ± + GG 


jets + x?Ai -> jets + 77 + GG 




leptons. The NLSP lifetime is a free parameter in GGM SUSY. Only prompt neutralino decays 
are considered in this analysis. 

Previous searches for gauge-mediated SUSY breaking were performed by the ATLAS exper- 
iment with 36pb _1 [151, l.lfb -1 El, and 4.8 ftT 1 02 of pp collision data, by CMS with 
36 pb _1 HU, as well as by experiments at the Tevatron (HUHl, LEP EHSU, and HERA (23. 

2.2 Simplified Models 

The experimental results of the single-photon and diphoton analyses are in addition interpreted 
in the context of Simplified Models (SMS) ||261 - |3T| . In SMS, a limited set of hypothetical particles 
and decay chains are introduced to produce a given topological signature such as the single or 
diphoton final state studied in this analysis. The amplitudes describing the production and 
decays of these particles are parametrized in terms of the particle masses. In particular, pairs 
of gluinos are initially produced that decay to jets and either a neutralino, and chargino or 
two neutralinos as shown in Fig. [2] The neutralino is then forced to decay into a photon and 
undetected LSP while the chargino is forced to produce a W boson resulting in either a single- 
photon or a diphoton final state. Simplified Models provide a benchmark different from other 
constrained models such as the GGM SUSY scenario for comparing different search strategies 
on a topological level. They also facilitate limit comparisons with other final state topologies. 

2.3 Universal Extra Dimensions 

Diphoton final states with large E™ 1SS similar to those expected from GGM SUSY scenarios are 
also predicted by models based on universal extra dimensions. Here the existence of additional 
compactified dimensions are predicted in which SM fields can propagate. The UED scenario 
provides several significant consequences including gauge-coupling unification, supersymme- 
try breaking, and other phenomena beyond those predicted by the standard model |3ll32|. The 
propagation of SM particles through the additional dimensions leads to the existence of a series 
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of excitations for each SM particle, known as a Kaluza-Klein (KK) tower, which can decay via 
cascades involving other KK particles until reaching the lightest Kaluza-Klein particle (LKP), 
which is the first level KK photon. SM particles such as quarks and leptons can also be pro- 
duced in the cascades. 

The UED space can be embedded in a larger space that has n large extra dimensions (LED) 
where only the graviton propagates with a (4 + n) -dimensional Planck scale (Mp) of a few TeV. 
In this case the LKP is allowed to decay gravitationally producing a photon and a graviton. As 
the dominant production mechanism at the LHC is from the strong interaction, KK quark and 
gluon pairs are produced, cascading down to two LKP decays resulting in the two photon plus 
jet(s) and E™ 1SS final state. Previous UED studies have been performed by the DO experiment 
at the Tevatron ||20| and most recently by ATLAS fl5|. 

3 The CMS Detector 

The central feature of the CMS detector is a superconducting solenoid, of 6 m internal diam- 
eter, providing an axial magnetic solenoid of 3.8 T along the beam direction. Within the field 
volume are a silicon pixel and strip tracker, a crystal electromagnetic calorimeter (ECAL), and 
a brass/ scintillator hadron calorimeter (HCAL). Charged particle trajectories are measured by 
the silicon pixel and strip tracker system, covering < (p < 2rc in azimuth and \rj\ < 2.5, 
where the pseudorapidity is tj = —In [tan 6/2], and 9 is the polar angle with respect to the 
counterclockwise-beam direction. Muons are measured in gas-ionization detectors embedded 
in the steel return yoke. Extensive forward calorimetry complements the coverage provided by 
the barrel and endcap detectors. 

The electromagnetic calorimeter, which surrounds the tracker volume, consists of 75 848 lead- 
tungstate crystals that provide coverage in pseudorapidity \tj\ < 1.479 in the barrel region (EB) 
and 1.479 < \rj\ < 3.0 in two endcap regions (EE). The EB modules are arranged in projective 
towers. A preshower detector consisting of two planes of silicon sensors interleaved with a 
total of 3 Xo of lead is located in front of the EE. In the region \rj\ < 1.74, the HCAL cells have 
widths of 0.087 in pseudorapidity and azimuth (</>). In the (rj, (p) plane, and for < 1.48, the 
HCAL cells map on to 5 x 5 ECAL crystal arrays to form calorimeter towers projecting radially 
outwards from close to the nominal interaction point. At larger values of \t]\, the size of the 
towers increases and the matching ECAL arrays contain fewer crystals. Within each tower, the 
energy deposits in ECAL and HCAL cells are summed to define the calorimeter tower energies, 
subsequently used to provide the energies and directions of hadronic jets. In the 2011 collision 
data, unconverted photons with energy greater than 30 GeV are measured within the barrel 
ECAL with a resolution of better than 1% 1133 1 . The HCAL, when combined with the ECAL, 
measures jets with a resolution AE/ E 100 %/ \JE [GeV] 5 %. The CMS detector is nearly 
hermetic, allowing for reliable measurements of E™ ss . A more detailed description of the CMS 
detector can be found in Ref. C34l . 

4 Data Selection 

The data sample used in this analysis was recorded during the 2011 pp run of the LHC at a 
center-of-mass energy of 7 TeV and corresponds to an integrated luminosity of 4.93 fb _1 . Events 
were selected using the CMS two-level trigger system requiring the presence of at least one 
high-energy photon and significant hadronic activity or at least two photons. The first level of 
the CMS trigger system, composed of custom hardware processors, uses information from the 
calorimeters and muon detectors to select events in less than 3.2 ^s. The High Level Trigger 



4.1 Photon and Electron Reconstruction and Identification 



5 



processor farm further decreases the event rate from around 100 kHz to around 300 Hz, before 
data storage. 

Photon triggers are utilized to select both the signal candidates and control samples used for 
background estimation. The efficiency for signal events to pass the trigger requirements ranges 
around 40-60%, while the efficiency for signal events which pass the photon offline selection is 
estimated to be greater than 99%. The single-photon search is based on the photon-Hx trigger 
requiring the presence of one photon with Ej > 70 GeV and the quantity Hj, the scalar sum 
of transverse momenta of reconstructed and calibrated calorimeter jets with pj > 40 GeV and 
\tj 1 < 3.0 in the event. Because of the continuous increase in the instantaneous luminosity, the 
trigger evolved with time from Hj > 200 to Hj > 400 GeV. An inefficiency of this trigger 
during a short time period of data taking restricts the single-photon analysis to an integrated 
luminosity of 4.62 fb _1 . The diphoton measurement using an integrated luminosity of 4.93 fb _1 
of pp collisions is based on a diphoton trigger with an Ej threshold of 36 GeV (22 GeV) for the 
leading (sub-leading) photon. 

4.1 Photon and Electron Reconstruction and Identification 

Photon candidates are reconstructed from clusters of energy in the ECAL. The photon iden- 
tification requires the ECAL cluster shape to be consistent with that expected from a photon, 
and the hadronic energy detected in the HCAL behind the photon shower not to exceed 5% 
of the ECAL energy. To suppress hadronic jets being misreconstructed as photons, we require 
photon candidates to be isolated from other activity in the tracker, ECAL and HCAL. A cone 
of AR = \J (At/) 2 + (A(p) 2 = 0.3 is constructed around the direction of the photon candidate, 
and the scalar sums of transverse energies of tracks and calorimeter deposits within this AR 
cone are determined, after excluding the contribution from the photon candidate itself. These 
isolation sums for the tracker, ECAL and HCAL are added to form J C omb- This combined isola- 
tion sum is corrected for contributions from additional pp interactions (pileup) II35II other than 
the hard scattering that produced the photon(s) and jets of interest I^°^ b = 4omb ~~ ^pileup- The 
corrected combined isolation is required to be ^°^ b < 6 GeV, which is based on an optimiza- 
tion of s/Vb as a figure of merit, where the signal S is from simulated SUSY-GGM events (see 
Section[5]) and the background B corresponds to a multijet simulated sample. As a cross check, 
data from a multijet-enriched sample consisting of events with low missing transverse energy 
£miss < 3Q Q e v / where the photon candidates passing all analysis requirements except the iso- 
lation cut, were also used as background sample. Using the same signal GGM sample, this test 
also results in an optimal value of ^°mb < ^ GeV. 

With increasing instantaneous luminosity during the 2011 LHC operation, the number of inter- 
actions per bunch crossing has also increased, resulting in an approximately linear rise in the 
occupancy of the detector. The energy in the photon isolation cone is sensitive to pileup effects. 
In an effort to reduce the dependence on the variation of pileup, an effective energy propor- 
tional to the amount of pileup Epii eu p = P x ^eff is subtracted from the combined photon- 
isolation variable. The p variable describes the amount of transverse momentum added to the 
event per unit area, e.g. by minimum bias particles, and A e a corresponds to an effective area 
determined from the slope of the average isolation energy versus p. p is described in more de- 
tail in Ref. 1351 , and it and the p x A e ff isolation compensation factor is calculated from the data 
on an event by event basis. Separate effective areas are calculated for the ECAL and HCAL 
isolation. 

The criteria above are efficient for the selection not only of photons but also of electrons. To 
reliably separate them, we search for hit patterns in the pixel detector consistent with at least 
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a single pixel hit matching a track from an electron. The candidates without pixel match are 
considered to be photons. Otherwise they are considered to be electrons, which are used to 
select control samples for background estimation. 

Photons that fail either the shower shape or combined isolation requirement are referred to 
as misidentified photons. These objects are predominantly electromagnetically-fluctuated jets 
and are used for the background estimation based on data. The definition of the misidentified 
photon is designed to be orthogonal to our real candidate photons, but still similar to that of 
the real photon definition to provide an accurate background estimate. An upper bound on 
^comb * s introduced in order to avoid events with highly isolated misidentified photon objects 
where the resolution on E™ 1SS is expected to be different compared to events with real photons. 
An upper cut of J™mb < 30 GeV (20 GeV) was found optimal for the single-photon (diphoton) 
analysis. 

Photons which convert in the tracker material ahead of the ECAL are reconstructed and counted 
as photon objects. These photons can have slightly higher isolation sums that unconverted pho- 
tons or, if they convert in the pixel detector, can be counted as electrons. Both possibilities of 
contamination have been studied and found to be negligible in this analysis. 

4.2 Jet and Missing Transverse Energy Reconstruction and Identification 

Jets and E™ 1SS are reconstructed with a particle-flow (PF) technique p6ll37| . The PF event re- 
construction consists of identifying every particle with an optimized combination of all sub- 
detector information. The energy of photons is obtained directly from the ECAL measurement, 
corrected for detector effects. The energy of electrons is determined from a combination of the 
track momentum at the primary interaction vertex, the corresponding ECAL cluster energy, 
and the energy sum of all bremsstrahlung photons attached to the track. The energy of muons 
is obtained from the corresponding track momentum. The energy of charged hadrons is deter- 
mined from a combination of the track momentum and the corresponding ECAL and HCAL 
energy, corrected for detector effects, and calibrated for the non-linear response of the calorime- 
ters. Finally, the energy of neutral hadrons is obtained from the corresponding calibrated ECAL 
and HCAL energy. 

All these particles are clustered into jets using the anti-Zcj clustering algorithm [38| with a dis- 
tance parameter of 0.5. The jet momentum is determined as the vectorial sum of all particle 
momenta in this jet and is found in the simulation to be within 5% to 10% of the true momen- 
tum over the whole pj spectrum and detector acceptance. An offset correction is applied to 
take into account the extra energy clustered in jets due to multiple pp interactions within the 
same bunch crossing, thereby reducing the dependence of jet energies on the instantaneous 
luminosity. Jet energy corrections are derived from simulation studies and are compared with 
in situ measurements using the energy balance of dijet and photon+jet events. Additional se- 
lection criteria are applied to each event. For example, jets identified to originate in spurious 
jet-like features from isolated electronic noise patterns in HCAL and ECAL are removed from 
the sample ||37| . 

Jets selected for this analysis are required to have transverse momentum pj > 30 GeV, \rj \ < 2.6 
and to satisfy the following jet-selection requirements. The neutral-hadron fraction as well as 
the electromagnetic fraction of energy contributing to the shower created by the jet should 
each be <0.99, and the charged hadron fraction is required to be greater than zero. Events 
must contain at least one jet isolated from the photon candidates by AR > 0.5 for the events to 
be retained in the signal sample. 
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Table 2: Summary of the signal and control sample selection criteria used for the single-photon 
and diphoton analyses. Electron (ee) and misidentified photon (ff) categories are used in back- 
ground estimations described in Sections [7] and [8] The exclusive bins in E miss are used in the 
limit setting procedure. 





Signal 


Single phot( 
Multijet control 


)n 

EWK control 


Signal 


Diphol 
ee control 


ton 

// control 


'comb [GeV] 
pixel seed 


< 6 
veto 


> 6, < 30 
veto 


< 6 
required 


< 6 
veto 


< 6 
required 


> 6, < 20 
veto 


Trigger 

Photon(s) 

PFJet(s) 

Hj 

£miss 


y-Hj trigger with 

p\ > 70 GeV, Hj > 400 GeV 

(using p J T ts > 40 GeV, \rj\ < 3.0) 

p\ > 80 GeV, \tj\ < 1.4 
p J T ts1 ' 2 > 30 GeV, \tj\ <2.6 
Hj > 450 GeV 

(using p^ 7 > 40 GeV, \tj\ < 3.0) 
E miss > 100 GeV (6 excl. bins in E miss ) 


77 trigger with 

pf' 2 > 36 (22) GeV 

Pj 1 ' 2 > 40 (25) GeV, \tj\ < 1.4 
p ] j > 30 GeV, \tj\ <2.6 

E miss > 50 GeV (5 excl. bins in E miss ) 



4.3 Single-Photon and Diphoton Event Selections 

The single-photon analysis requires the scalar sum of the transverse energy of all jets and all 
photons in the event, Hj, to be larger than 450 GeV, where the photon- Hj trigger is fully ef- 
ficient. To closely resemble the trigger requirement, calorimeter jets with pj > 40 GeV and 
I tj\ < 3.0 are used for the Hj calculation, but with the addition that these jets are pileup cor- 
rected. Both real and misidentified photons are included in the Hj calculation. Since the photon 
objects are also reconstructed as jets, the pj of the jet is used in the Hj calculation instead of the 
photon object, if the transverse momentum ratio between jet and photon object is greater than 
95% and the photon and jet are within AR < 0.3. This avoids a bias in Hj and E miss due to the 
different isolation requirements for the genuine photon candidates and the misidentified pho- 
tons in the multijet control samples. In addition, a photon with Ej > 80 GeV within \tj\ < 1.4 
and at least two jets with pj > 30 GeV and \rj\ < 2.6 are required. Events with isolated leptons 
are not rejected, and the lepton momenta are not included in the Hj determination. 

To be within the full efficiency of the 77 trigger with an Ej threshold of 36 GeV (22 GeV) on 
the leading (sub-leading) photon, the diphoton offline analysis requires at least two photons 
with Ej > 40 GeV (25 GeV) for the leading (sub-leading) photon in the event and at least one 
jet with pj > 30 GeV and \t]\ < 2.6. Table [2] contains a summary of the signal sample selection 
criteria for the single-photon and diphoton analyses. It also includes information on the back- 
ground control samples described in Section [6] as well as the search region for new physics in 
the variable of transverse missing energy as further discussed in Section [9] 

5 Simulated Samples 

Although this analysis uses methods based on data to estimate the background components, 
simulated samples are used to evaluate less significant backgrounds, which might be difficult 
to measure directly from the data, or to model the new physics (NP) signals and to validate the 
performance of the background estimation from data. 
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5 Simulated Samples 



Table 3: Parameters varied in GGM signal scans used in the interpretation. Grid values along 
either axis in the scan are offset by 10-20 GeV to prevent degeneracies between the generated 
particles. 



Scan name 


Squark mass 


Gluino mass 


Bino mass 


Wino mass 


Squark-Gluino (Bino) 
Squark-Gluino (Wino) 
Gluino-Bino 
Gluino-Wino 
Wino-Bino 


400-2000 GeV 
400-2000 GeV 
5000 GeV 
5000 GeV 
5000 GeV 


400-2000 GeV 
400-2000 GeV 
300-1500 GeV 
300-1000 GeV 
5000 GeV 


375 GeV 
5000 GeV 

50-1500 GeV 
5000 GeV 

5-1000 GeV 


2000 GeV 
375 GeV 
2000 GeV 
100-1000 GeV 
115-1000 GeV 



The simulated samples used in this search are produced in several ways. Depending on the 
process either the pythia (39| or MadGraph [40] Monte Carlo (MC) event generators are 
used to generate event kinematics and fragment partons into jets. For most simulated data, in 
particular to study SM backgrounds, the generated events are passed through the full Geant4- 
based [41 j CMS detector simulation. Because of the large number of individual simulated sam- 
ples required in the NP parameter space scans used in the interpretation of results in the light 
of NP, a fast detector simulation [42J based on a full description of the CMS detector geometry 
and a parameterization of single-particle showers and response is utilized to reduce the com- 
putation time for those samples. Event pileup corresponding to the luminosity profile of the 
analyzed data is added to all simulated samples and the generated events are reconstructed 
using the same software program as for the collision data. 

In interpreting our results, multiple samples of simulated signal data are produced by varying 
model parameters individually (as in the case of the UED interpretation) or in pairs (in the case 
of the GGM and SMS interpretations). General gauge-mediated SUSY breaking requires the 
LSP to be a gravitino, and the NLSP needs to be a wino-like or bino-like neutralino to produce 
a final state with photon(s) plus E™ 1SS . Bino-like neutralinos decay most of the time into a pho- 
ton. Wino-like neutralinos decay mostly into Z bosons, but they also decay into a photon ~20% 
of the time, allowing our measurement to be sensitive to this channel. In the GGM scans, other 
SUSY particles are decoupled (forced to have high mass) in order to leave only the possibility 
of light squarks, gluinos and the desired neutralino NLSP or neutralino /chargino co-NLSP as 
kinematically allowed production particles. Table [3] shows the mass parameters varied in the 
five GGM planes investigated in this analysis |fl3|. The masses of these particles take values 
within the ranges indicated in the table as different scan grids are produced. In particular, 
the SUSY mass spectra are calculated in form of files following the SUSY Les Houches Ac- 
cord (SLHA) [43] utilizing SuSpect El with decay tables from SDECAY [45 1. The SUSY GGM 
events are generated in a three-dimensional grid of squark, gluino, and NLSP masses. Squarks 
are taken to be degenerate in mass and all other SUSY particles are assumed to be heavy. In the 
scans where the NLSP mass is varied, the "next-to-next to LSP" (usually a gluino) is required 
to have a higher mass, resulting in scans that only span above the diagonal in the correspond- 
ing mass plane. This is also the case for the Simplified Model scans described below. In the 
"Wino-Bino" scan shown at the bottom of Table |3j we decouple the squarks and gluinos, leav- 
ing only electroweak production of wino-like neutralino/ charginos to study our sensitivity to 
electroweak production of SUSY. 

For the Simplified Model interpretation, more controls are exerted over the production and 
decay of sparticles, which are often forced to decay into a certain final state, e.g., 100% of the 
time. Two parameter scans referred to as the W7 SMS (Fig. |2] left) and the 77 SMS (Fig. |2| 
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right) are used in this analysis. They both span a grid in gluino and neutralino/ chargino mass 
space, forcing the initial pair production of gluinos, which then decay to jets and neutralino or 
chargino. In the 77 Simplified Model, both gluinos are forced to decay to jets and neutralinos, 
which in turn decay to photons. The W7 SMS forces one gluino to decay to a chargino, which is 
forced to always produce a W boson, and the other gluino decays as in the 77 Simplified Model. 
The 77 scan produces final states to which both the single-photon and diphoton analyses are 
sensitive, while the W7 SMS scan is interpreted only through the single-photon analysis. The 
production cross sections of the GGM and SMS scans [46J are calculated at next-to-leading 
order (NLO) plus next-to-leading log in QCD using the PROSPINO program HZrEH- Except 
for the GGM "Wino-Bino" scan, the production in these scans is dominated by gluino-gluino, 
gluino-squark, and squark-squark production. 

Simulated signal samples for the UED interpretation are generated using the UED model as 
implemented at leading order (LO) in PYTHIA [39 j. Parameters for the UED model investigated 
in this analysis including the LO cross section are chosen to match previous UED searches by 
other experiments |[T5ll20|| . The UED model has two varying parameters, the ultraviolet cutoff 
A and the radius of compactification R. In this study R is chosen as a free parameter while 
A is set to satisfy the relation AR = 20 £Q. Additional parameters that are used in the MC 
generation of the signal are chosen as follows. The number of large extra dimensions is N = 2 
or 6, the (N + 4) -dimensional Planck scale Md is 5TeV, while the number of KK excitation 
quark flavors is five. Sample points of 1/R ranging from 900 to 1600 GeV are produced in 
increments of 50 GeV. 



6 Background Estimation Methodology 

The NP signature of the photon(s) plus E™ 1SS final state can be mimicked by SM processes 
in several ways. The largest backgrounds are due to events without true E™ ss resulting from 
abundant hadronic processes, such as direct photon plus jets processes, and multijet production 
with electromagnetically rich jets misidentified as photons, which result in events with the 
same topology as the NP signal. The missing Ej in these hadronic events comes from poorly 
measured hadronic activity in the event. This background is referred to as background with 
false E™ 1SS or as QCD background. The E™ 1SS resolution for this background is much poorer 
than the resolution of the total Ej of the photon(s) and is determined by the resolution of the 
hadronic energy in the event. The strategy for determining the shape of the E™ 1SS distribution 
for the QCD background is to find a control sample that reproduces the hadronic activity in 
the candidate sample while having no significant true E™ 1SS that mimic a substantial missing 
Ej contribution. 

The second kind of background comes from processes with true Eij? iss . It is dominated by W7 
events and W plus jets production where the W decays into an electron plus a neutrino, with 
the electron or jet misidentified as a photon and the neutrino leading to E™ 1SS . We refer to this 
sample as background with true E™ 1SS or electroweak (EWK) background and it is determined 
in the following way. Since the photon is expected to behave almost identically to an electron 
in the electromagnetic calorimeter, electrons can be mistaken as photons except that electrons 
have hits matching the particle track in the pixel detector. We measure the electron-photon 
misidentification rate / e ^7 and determine the con tribution of the EWK background by apply- 



ing /e^7 to our EijP lss distribution (see Section 6.1 1. The rates of other processes with true Ej 



that have single photon or diphotons in their final states are quite small and are discussed for 
the single-photon and diphoton analyses separately in Sections [7] and [H] 
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7 Single-Photon Analysis 



6.1 Electron Misidentification Rate 

We determine the probability to misidentify an electron as a photon, by fitting the mass of 
the Z — > e + e~ peak seen in the ee and e7 mass spectra, and comparing the integrals of these 
fits. For this purpose we identify a sample of ee events where pixel matches are required on 
both objects that otherwise satisfy the photon selection requirements (see details of diphoton 
analysis in Section [8]). The sample has the same requirements imposed on it as the real 77 
sample, except a pixel match is required for one of the electromagnetic objects. 

We extract the electron misidentification fraction from the ee and e7 spectrum using the num- 
ber of observed Z — > ee events in the ee mass spectrum given as N ee = (1 — / e ^ 7 ) 2 Nz true 
where Nz true is the number of true Z — > ee events. The observed Z — > ee peak in the e7 spec- 
trum is N e7 = 2 \f e ^ 7 (l — / e ^ 7 )] Nz true leading to / e ^ 7 = N e7 / (2N ee + N e7 ). We can calculate 
the number of Z — > ee events expected in the 77 spectrum using N 77 = (/ e ^ 7 ) 2 x N ee / (1 — 
fe^y) 2 and cross check the number of observed diphoton events. 

We measure / e ^ 7 in bins of photon transverse momentum. The overall misidentification rate 
integrated over the whole p T range is determined as / e ^ 7 = 0.015 ± 0.002 (stat.) ± 0.005 (syst). 
This number is used for the diphoton analysis, while for pj > 80 GeV a misidentification rate 
of / e -^ 7 = 0.0080 ± 0.0025 (stat.) is determined. The latter rate is used for the single-photon 
search since pj(j) > 80 GeV is the momentum region relevant for this analysis. 

7 Single-Photon Analysis 

The single-photon analysis targets especially SUSY scenarios in which the lightest gaugino 
comprises a large non-bino-like mixture. In this case the branching fraction of the lightest 
gaugino to a photon and the gravitino LSP is reduced and decays into other bosons like W, Z, 
or Higgs occur, leading to additional jets and possibly leptons in the final state, suppressing 
events with more than one photon. Events with leptons or more than one photon are not 
removed in the single-photon analysis. The potential overlap with the diphoton selection has 
been studied and is found to be negligible. 

7.1 Background Estimation 

The dominant background in the single-photon analysis is a composition of processes such as 
7+jets and multijet QCD production with one jet misidentified as a photon. The shape of the 
E™ 1SS distribution is similar for both background contributions, as the event topologies are very 
similar. Therefore, these two contributions to the QCD background are estimated together from 
the same control sample. This background sample is selected by applying the signal selection 
requirements, except that the photon candidate is required to fail the photon identification cri- 
teria but to satisfy a loose isolation requirement. Such misidentified photon candidates follow 
a definition orthogonal to the photon identification criteria in the signal selection. The back- 
ground control sample is weighted to correct for the difference in pj spectra of misidentified 
and genuine photons. The weights as a function of the photon transverse energy are deter- 
mined in bins of pj from the ratio of events in the misidentified and genuine photon samples 
for E™ 1SS < 100 GeV, which is taken as a signal-depleted region for the normalization of the 
QCD background to the single-photon data. 

The EWK background contribution is much smaller than the QCD background. The dominant 
contributions are from tt production or events with W or Z bosons with one or more neutrinos 
in the final state in which the electron is misidentified as a photon. This background is modeled 
from the data using an electron control sample selected by the same trigger as the signal dataset. 



7.2 Results 
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Figure 3: Missing Ej spectrum of single-photon data (dots with error bars) compared to various 
SM background predictions (solid colored histograms). The shaded area indicates the uncer- 
tainty in the total background prediction. The E™ 1SS spectrum for two example GGM points 
(red upper and blue lower solid curves with masses of m^/m^/m^ inGeV) on either side of 
our exclusion boundary are also shown. At the bottom, the ratio of data over standard model 
prediction is shown as a function of E™ 1SS . The error bars take into account only the statistical 
error of the data sample, while the hatched area is the uncertainty in the expected background 
from the SM processes. 



The electron control sample is weighted according to the misidentification rate, f e ^y, measured 



in Z — > ee events, as discussed in Section 6.1 



Additional backgrounds can contribute due to initial-state radiation (ISR) and final-state radi- 
ation (FSR) of photons. Both ISR and FSR, in events with electrons in the final state, are already 
covered by the EWK background prediction from data. The remaining contributions from W, 
Z, and tt events are taken from MC simulation. 

7.2 Results 

The dominant systematic uncertainty in the background estimation arises from the small num- 
ber of events in the misidentified-photon control sample. The statistical uncertainty associated 
with the E™ 1SS < 100 GeV sample, where the normalization of misidentified and genuine pho- 
ton samples is calculated in bins of photon pj, is propagated as a systematic uncertainty. The 
uncertainty is taken to be correlated among the E™ 1SS bins, as a given E™ 1SS bin receives contri- 
butions from several photon pj normalization bins. The method assumes, that the E™ 1SS and the 
photon momentum are uncorrelated. This has been validated in simulation up to 5%, which is 
assigned as additional systematic uncertainty. 

In comparison, the systematic uncertainty due to the statistically limited electron control sam- 
ple used for the electroweak background prediction is negligible. In addition, the small un- 
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8 Diphoton Analysis 



Table 4: Resulting event yields for the >1 photon and >2 jet selection in 4.62 fb 1 of data for 
six distinct signal search bins. 



E™ 1SS bins [GeV] 


100-120 


120-160 


160-200 


200-270 


270-350 


> 


350 


QCD (from data) 


262 ±37 


173 ±27 


82 ±24 


55 ±14 


29 ±11 


6.8 


±4.2 


e— > 7 (from data) 


4.5 ±1.9 


6.0 ±2.5 


3.2 ±1.3 


2.3 ±1.0 


0.8 ±0.4 


0.4 


±0.2 


FSR/ISR (W,Z) 


4.7 ±1.3 


8.2 ±1.8 


5.5 ±1.5 


5.4 ±1.3 


4.0 ±1.3 


1.7 


±0.9 


FSR/ISR (tt) 


0.6 ±0.3 


1.7 ±0.6 


0.9 ±0.4 


0.5 ±0.4 


0.4 ±0.3 


< 


0.01 


Total SM estimation 


272 ±37 


189 ±27 


91 ±24 


63 ±14 


34 ±11 


8.8 


±4.3 


Data 


283 


199 


70 


39 


20 
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certainty in the electron misidentification rate f e ^j = 0.008 ± 0.0025 is propagated resulting 
in small systematic uncertainties in the EWK background prediction. Finally a conservative 
uncertainty of 50% on the ISR and FSR contributions to the W/ Z and tt cross sections is added. 

All background components are shown in Fig. [3] together with the data (points with errors 
bars) and two GGM benchmark signal samples, one excluded (red line) and one not excluded 
(violet line) by this analysis. The same information is summarized in Table[4] No excess beyond 
standard model predictions is observed. 

8 Diphoton Analysis 

The diphoton analysis is most sensitive to SUSY scenarios in which the lightest neutralino is 
bino-like decaying into a photon and the gravitino as LSP, as well as models predicting uni- 
versal extra dimensions. To keep the analysis as inclusive as possible, no veto is applied on 
additional leptons in the event. 

8.1 Background Estimation 

To estimate the QCD background from data in the diphoton analysis, two different datasets are 
utilized. The first sample contains two misidentified photons, and in what follows referred to 
as the // ("fake-fake") sample, comprising multijet events. This is the main dataset to estimate 
the QCD background. The second data sample contains events with two electrons (ee) with an 
invariant mass between 81 and 101 GeV, and is dominated by Z — > ee decays. The ee sample is 
used to study systematic effects on our background estimate. 

Comparing the E™ 1SS resolution between diphoton signal and background events, the Ej res- 
olution for electrons and misidentified photons is similar to the resolution for true photons. 
It is negligible compared with the resolution for the hadronic energy, which dominates the 
£miss resolution. The events in both control samples are reweighted to reproduce the diphoton 
transverse energy distribution in the signal data sample, and, therefore, the transverse energy 
of hadronic recoil against the photons. The E™ 1SS distributions in the reweighted control sam- 
ples show good agreement with the diphoton signal samples within uncertainties as shown 
for the // sample in Fig. |ij The shape of the E™ 1SS distribution for the // sample is used to 
determine the magnitude of the QCD background after normalizing the // background shape 
to the diphoton data in the region of low missing transverse energy E^ 1SS < 20 GeV, which is 
dominated by QCD background. We choose to use the prediction from the // sample as the es- 
timator of the QCD background while the difference from the sideband-subtracted ee sample to 
the // estimate is taken as an estimate of the systematic uncertainty in the determination of the 
QCD background. The ee sample has been corrected for a small contribution from diboson pro- 
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Table 5: Number of diphoton candidates from data as well as estimates of QCD and EWK 
background in bins of E™ ss . The first error is statistical and the second is systematic for each 
entry. 



EJf 1SS bins [GeV] 


50-60 


60-70 


70-80 


80-100 


> 100 


QCD background 


183.8 ± 17.7 ± 12.5 


67.3 ±10.7 ±13.6 


15.4 ±5.1 ±11.5 


9.4 ±4.0 ±0.7 


10.1 ±4.2 ±1.4 


EWK background 


6.5 ±0.3 ±2.2 


3.1 ±0.2 ±1.0 


2.2 ±0.2 ±0.7 


2.2 ±0.2 ±0.8 


2.9 ±0.2 ±1.0 


Total background 


190.3 ±17.7 ±12.7 


70.4 ±10.7 ±13.7 


17.6 ±5.1 ±11.5 


11.6 ±4.0 ±1.0 


13.0 ±4.2 ±1.7 


Data 


199 


63 


26 


26 
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duction (WZ and ZZ) using PYTHIA with NLO cross section resulting in a correction of 0.2-18% 
depending on E^ 1SS bins, and ee events with true EJp lss . As an illustration of the reliability of the 
QCD background estimate, in the E™ 1SS control region from 30 to 50 GeV, 3443 candidate dipho- 
ton events are observed in the sample requiring >1 jet in the event. In the same E™ 1SS region 
the prediction from the // and ee sample yields 3636 ± 79 (stat.) and 3045 ± 26 (stat.) events, 
respectively. 

The estimated EWK background is determined with the ee and samples as described in 
Section|6]and is calculated to be much smaller than the QCD background. Other backgrounds 
such as Z77 — > 1/1/77, W77 — > £1/77, tt77, or Z77 events where the Z — > tt is followed by a r 
decay such as T — > nv or t — > e{}i)vv have been found to be <0.1% using simulations. 

Drell-Yan events can also contribute as background if both electrons are misidentified as pho- 
tons. While the Drell-Yan process does not have true E™ 1SS , it can have mismeasured E™ 1SS 
due to resolution effects in the accompanying hadronic activity. Given the high expected elec- 
tron pixel match efficiency, and the relatively low cross section for Drell-Yan production, the 
contribution from this background is also negligible. 

8.2 Results 

The E^ ss distribution in the 77 sample requiring > 1 jet in the event is presented in Fig. E] 
as points with errors bars. The green shaded area shows the estimated amount of the EWK 
background while the QCD background prediction from the // sample is shown in grey after 
normalization to the 77 sample minus the estimated EWK contribution in the region E™ ss < 
20 GeV. The hatched areas indicate the total background uncertainties. 

Table § summarizes the observed number of 77 events in bins of E™ 1SS as well as the expected 
QCD and EWK background with statistical and systematic uncertainty. The systematic error is 
determined from the difference between the // sample used to predict the QCD background 
and the ee sample utilized as an alternative background estimate after the ee data are sideband 
subtracted and corrected for a small diboson contributions. For the region of large missing 
transverse energy, no excess of data over the SM expectation is found. We observe 11 diphoton 
events with E™ ss > 100 GeV while the total background expectation is calculated to be 13.0 ± 
4.2 (stat.) ± 1.7(syst.) events. 

9 Interpretation in Models of New Physics 

We determine the efficiency for NP signal events to pass our analysis selections by applying 
correction factors derived from data to the MC simulation of the signal. Since there is no large 
clean sample of genuine photons in the data, we rely on the similarities between the detector 
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Figure 4: The E™ 1SS spectrum of 77 data compared to QCD prediction together with small EWK 
background for events with at least one jet. The hatched areas indicate the total background 
uncertainties. Two example GGM points (dashed red upper and solid blue lower curves with 
masses of m^/ m^/ m^o in GeV) on either side of our exclusion boundary are also shown. At the 



bottom, the ratio of data over standard model prediction is shown as a function of E2p iss . The 
error bars take into account only the statistical error of the data sample, while the hatched area 
is the error on the expected background from the SM processes. 
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response to electrons and photons to extract the photon identification efficiency. A scale factor 
is obtained and applied to the photon efficiency in MC simulation by forming a ratio between 
the electron efficiency from Z — > ee events that pass all photon selections (except for the pixel 
match) and the corresponding electron efficiencies from simulation. The obtained data-to-MC 
scale factor 0.994 ± 0.002 (stat.) ± 0.035 (syst.) is applied to the photon efficiencies obtained from 
MC simulation. Other sources of the larger systematic uncertainties in the signal yield include 
the error on integrated luminosity (2.2%) Il52l , pileup effects on photon identification (2.5%), 
and small parton distribution functions (PDF) uncertainties in the acceptance. Systematic un- 
certainties in the theoretical cross section prediction consist of the PDF uncertainty (4-66%) and 
renormalization scale (4-28%) uncertainty depending on the parameters of the NP signal. 

The goal of this analysis is to find evidence for the production of NP by observing an excess 
of events above the SM background in the high E™ 1SS region of the single-photon and diphoton 
signal samples. Since no such excess is observed, upper limits are derived on potential signals 
of various NP models. The statistical approach used to derive limits constructs a test statistic 
as the product of likelihood ratios in bins of E™ 1SS . These likelihoods are functions of the pre- 
dicted signal and background yields in each bin. Systematic uncertainties are introduced as 
nuisance parameters in the signal and background models. Log-normal distributions are taken 
as a suitable choice for the probability density distributions of the nuisance parameters in or- 
der to incorporate uncertainties in the background rates, integrated luminosity, and the signal 
acceptance times efficiency. 

In order to compare the compatibility of the observed data with a NP signal hypothesis, we 
use a LHC-style profiled likelihood test statistics [53 j. In particular, for the comparison of the 
data to a signal-plus-background hypothesis, where the signal and background expectations 
are functions of nuisance parameters 6 and the signal is scaled by a signal strength param- 
eter ]i, we construct a one-sided test statistic —Ihcicfp based on the profile likelihood ratio 
cfjt = C (data | }i,Q-^)l C (data \%v) with constraint < p < \i (53l. Here, 9 }l refers to the con- 
ditional maximum likelihood estimators of 6 given the signal strength parameter ji and the 
actual data. The pair of parameter estimators 9 and ft correspond to the global maximum of the 
likelihood. The modified frequentist CL S criterion Il54ll55l is used to determine upper limits on 
the cross section of a possible NP signal at the 95% confidence level (CL). 

To achieve optimal sensitivity, the limits are calculated in distinct bins and multiple exclusive 
search channels in E™ 1SS are combined into one test statistic considering the bin-to-bin cor- 
relations of the systematic uncertainties. For the single-photon analysis, six distinct bins for 
E™ iss > 100 GeV are used, [100,120), [120,160), [160,200), [200,270), [270,350), and [350,oo) given 
inGeV, while the diphoton analysis uses the following E™ ss ranges given in GeV: [50,60), 
[60,70), [70,80), [80,100), and [100,oo). These bins in E™ ss correspond to the event yields given 
in Tables [i] and [5] In general, the sensitivity is dominated by the highest E™ ss bin. Since both 
searches slightly overestimate the observed data in the highest E™ 1SS bin, the observed limits are 
generally slightly stronger than the expected limits. Some regions of the possible signal phase 
space, e.g. where the LSP receives only a small amount of transverse momentum, resulting in 
small E™ 1SS , also benefit from other search bins and therefore from the combination. 

A possible contamination by signal in the control samples used for the background estima- 
tion has been studied and was found to be negligible for the diphoton final state. For the 
single-photon analysis the expected contamination for a given signal is considered in the limit 
calculation in the signal-plus-background hypothesis. The background overestimation due to 
the contamination is typically a few percent, if the signal cross section is of the same order than 
the cross section limits. 
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9 Interpretation in Models of New Physics 



9.1 General Gauge Mediated SUSY Breaking 

Since the physical neutralinos x° and charginos x are an admixture of gaugino eigenstates, 
different scenarios of gaugino mixing have been studied. In the first case, referred to as bino- 
like, the lightest neutralino is assumed to be pure bino-like, while the lightest chargino is as- 
sumed to be heavy and decoupled. In this case, the production of the neutralino occurs mostly 
in the cascade decays of the squarks and gluinos, since the neutralino pair production cross 
section is very small. In the second case, referred to as wino-like, the neutralino and chargino 
have comparable mass and are assumed to be pure wino-like. In this case, both the neutralino 
and chargino are produced in squark and gluino decays, but direct chargino-neutralino pro- 
duction may also contribute. Furthermore, in the wino-like case, the expected event yields 
for the single-photon and diphoton analyses are reduced since the chargino (neutralino) may 
decay to a W (Z) and the gravitino (see Fig-[l). 

The resulting upper limits on the GGM production cross section, at 95% CL, as well as exclusion 
contours are shown in Fig. [5] for the gluino versus squark mass plane from 400 to 2000 GeV in 
squark and gluino mass, with the neutralino mass set at 375 GeV. This mass value is chosen 
to represent a reasonably light NLSP, but high enough to be outside current exclusion limits. 
For the wino-like scenario, the single-photon cross section upper limit is of order 0.003-0.1 pb at 
95% CL with a typical acceptance of ~7%. For the bino-like scenario, the diphoton cross section 
limit is of order 0.003-0.01 pb at 95% CL with a typical acceptance of ~30% for E™ ss > 100 GeV. 
Squark and gluino masses up to about 800 GeV are excluded in the wino-like scenario by the 
single-photon search, while the diphoton analysis excludes squarks and gluinos up to masses 
of ~1 TeV for a bino-like neutralino, both limits at 95% CL The corresponding 95% CL limits on 
the signal cross section and exclusion contours for the single-photon (diphoton) analysis in the 
bino-like (wino-like) scenario are available in Appendix |A| 

As further interpretation of the single-photon and diphoton results, Fig. [6] shows the exclusion 
contours in the plane of gluino versus neutralino mass for the single-photon wino-like and 
the diphoton bino-like scenarios. The diphoton search excludes gluino production for a bino- 
like neutralino up to about 1 TeV rather independent of the neutralino mass. The 95% CL upper 
limits on the signal cross section for the single-photon (diphoton) wino-like (bino-like) scenario 
in the gluino-neutralino mass plane as well as the corresponding single-photon bino-like and 
diphoton wino-like 95% CL limit plots and contours can be found in Appendix |A| 

Finally, we study for the first time in the final state with photons the electroweak production 
of winos, i.e. the pair and associated production of wino-like neutralinos and charginos, that 
decay to a bino-like NLSP by decoupling the squarks and gluinos leaving only electroweak 
production in the simulated samples. Figure [7] shows limits on the signal cross section and ex- 
clusion contours in the plane of wino-like versus bino-like gaugino mass for the single-photon 
and diphoton analyses, where the diphoton search excludes wino masses up to about 500 GeV 
almost independent of the bino mass. Since no continuous exclusion contour line can be drawn 
for the single-photon analysis, we can only present the 95% CL upper limits on the signal cross 
section. The corresponding 95% CL upper limits on the signal cross section in the plane of 
wino-like versus bino-like gaugino mass for the diphoton analysis are available in Appendix[A| 

9.2 Simplified Models 

In this section we interpret the results of our single-photon and diphoton search in terms of 
Simplified Models, which allow a presentation of our exclusion potential in the context of a 
larger variety of fundamental models, not necessarily in the GGM framework. For the SMS in- 
terpretation, we force the initial pair production of gluinos, which decay to jets and a neutralino 
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Figure 5: Observed upper limits at 95% CL on the signal cross section (left) and correspond- 
ing exclusion contours (right) in gluino-squark mass space for the single-photon search in the 
wino-like scenario (top) and the diphoton analysis for a bino-like neutralino (bottom). The 
shaded uncertainty bands around the expected exclusion contours correspond to experimental 
uncertainties, while the NLO renormalization and PDF uncertainties of the signal cross section 
are indicated by dotted lines around the observed limit contour. 



18 



9 Interpretation in Models of New Physics 



5. 1200. 

<D - 

£ 1100r 
1000^ 
900 r 
800- 
700 E- 
600 F 
500 
400 
300 



CMS Vs = 7 TeV 



>1y, >2 jets 



Wino-like x NLSP 

Observed 

Observed +1o theory 
Expected +1o exp. 



T 




57 2000 



— 1800 



^ 1600 
1400- 
1200 
1000 
800 
600 
400- 



CMS \!s = 7 TeV 



= 4.93fb , >2y's, > 1 jet 



I i I | I r— 
Bino-likex°NLSP 

Observed 

Observed ±1o theory 

Expected +1a exp. 




500 



1000 



M bin „ [GeV] 



Figure 6: Exclusion contours at 95% CL in the plane of gluino versus neutralino mass for the 
single-photon search in the wino-like scenario (left) and the diphoton analysis for a bino-like 
neutralino (right). 
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Figure 7: 95% CL exclusion contour and corresponding observed and expected contours in the 
bino-like versus wino-like gaugino mass for the diphoton (right) and the cross section limit for 
the single-photon analysis (left). 
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Figure 8: Results for Simplified Models in form of 95% CL upper limits on the cross section 
plus overlaid exclusion contours for the single-photon analysis in the W7 Simplified Model 
(left) and for the diphoton analysis in the 77 SMS interpretation (right). 

or chargino. Two cases are studied. Firstly in the 77 Simplified Model both gluinos decay to 
jets and neutralino, which are forced to decay to photons plus gravitino (see Fig. [2]) producing 
a final state with two photons. This model is sensitive to both the diphoton and single-photon 
analyses. Secondly in the W7 SMS, one gluino is forced to decay to a chargino, which always 
produces a W boson, and the other gluino decays as in the 77 SMS scan resulting in a photon, 
allowing only the single-photon analysis to be interpreted within this Simplified Model. 

The results in the form of upper limits on the cross section and overlaid exclusion contours, at 
95% CL, in the neutralino versus gluino mass plane are shown in Fig. [Ejfor the single-photon 
analysis in the case of the W7 Simplified Model, and for the diphoton analysis in the 77 SMS 
interpretation. The Simplified Model results in the gluino-neutralino mass plane are similar 
to the GGM interpretation resulting in slightly more stringent but similar limits as compared 
to the single-photon and diphoton contours shown in Fig. [6] This is not unexpected since 
both processes probe very similar production and decay chains and by construction, the SMS 
captures the main features of the full GGM model well. Additional figures such as the cor- 
responding acceptances in the gluino-neutralino mass plane for the single-photon (diphoton) 
analysis in the W7 (77) SMS interpretation and corresponding results from the single-photon 
analysis in the 77 Simplified Model are available in Appendix |A| 

9.3 Universal Extra Dimensions 

Diphoton final states with large £™ 1SS are also predicted by UED models [3| postulating the 
existence of additional spatial dimensions of compactification radius R. For the investigated 
model the UED space is embedded in an additional space of large extra dimensions where 
only the graviton propagates and the LKP decays gravitationally producing a photon and a 
graviton. The diphoton analysis results can thus be interpreted in the context of the UED 
model. The model parameters are chosen to match a study by the DO collaboration, which 
excludes 1 / R < A77 GeV II20II and a more recent result by the ATLAS experiment excluding 
1/R < 728 GeV [15|. To determine the effect of the number of large ED, n on the potential 
limit for UED, n was varied. By changing the number of large ED, the branching ratios of the 
different decay channels are changed but the overall UED production cross section remains the 
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Figure 9: Upper limit on the UED model cross section for n=6 at 95% CL compared with 
expected UED production cross sections (black diagonal line). The shaded region shows the 
uncertainty due to PDFs and renormalization scale on the expected limit. 



same. For n > 3 decays involving a heavy graviton with mass of order (1/R) dominate while 
for n = 2 decays involving light gravitons are more prevalent II56II . For n equal to 4 and 6, the 
£miss distributions are very similar allowing the comparison only for n = 6 to n = 2 where the 
£miss distribution is flatter resulting in a slightly lower efficiency. 

To determine the acceptance times efficiency UED signal simulated samples generated with 
1/R ranging from 900 to 1600 GeV as described in Section [5] are analyzed adopting the same 
selection criteria as used for the GGM diphoton analysis. The cross section upper limit for the 
production of KK particles, which would indicate the presence of UED, can be calculated in 
the same way as for the GGM limit calculation. The maximum UED production cross section 
is computed using the acceptance times efficiency from signal Monte Carlo simulations and 
the same luminosity, background estimate, and number of observed 77 signal events as for 
the GGM limit calculation. The signal acceptance times efficiency is rather flat in the region 
of interest ranging from about 0.42 at 1/R ~ 900 GeV to 0.46 at 1/R ~ 1600 GeV. The UED 
cross sections and the 95% CL upper limit on the signal cross section are interpolated and their 
intersection is determined and shown in Fig. |9j Uncertainties due to PDFs and renormalization 
scale are shown as the shaded region, while the intersection of the central value implies that 
the range of 1 /R < 1380 GeV for n = 6 is excluded with an expected limit of 1350 GeV. This is 
the best UED limit to date. For n = 2 the exclusion limit is reduced to 1350 GeV for an expected 
limit of 1340 GeV. The corresponding UED acceptance times efficiency distributions for n = 2 
and 6 as well as the 95% CL cross section upper limit for n = 2 are available in Appendix |A| 

10 Conclusions 

In summary, a search for physics beyond the standard model has been performed in single- 
photon and diphoton events using the E™ 1SS spectrum comparing data and SM background 
expectations. This search is based on 2011 CMS data comprising 4.93 fb -1 of pp collisions at 
^Js = 7 TeV. No evidence of NP is found and upper limits are derived for three theoretical 
interpretations. First, in the SUSY GGM model the single-photon (diphoton) analysis derives 
exclusion regions for the production cross section in the parameter space of squark and gluino 
masses of order 0.03 - 0.1 pb (0.003 - 0.01 pb) at the 95% CL for a wino-like (bino-like) scenario, 
corresponding to the exclusion of squark and gluino masses up to masses of order 800 GeV 
(1 TeV). Exclusion contours at the 95% CL are presented in the plane of gluino versus neu- 
tralino mass for a wino-like (bino-like) neutralino with the single-photon (diphoton) analysis. 
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In addition, for the first time, electroweak production is studied in the plane of wino-like versus 
bino-like gaugino mass where the diphoton search excludes wino masses up to ~ 500 GeV. 

The single-photon and diphoton analyses are in addition interpreted in the context of Simpli- 
fied Models resulting in similar exclusion limits and contours. Finally, the diphoton analysis is 
reinterpreted as a search for universal extra dimensions, leading to 95% exclusion values of the 
inverse compactification radius 1/R < 1380 GeV for n = 6 large extra dimensions constituting 
the currently best limit on the considered UED model. 
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A Supplemental Material 

The appendix contains additional figures such as limit contours from the three interpretations 
(GGM, SMS, and UED) that are not part of the main body of the paper but are submitted as 
supplemental material to the journal. 

A.1 GGM Interpretation 

This section contains additional 95% CL upper limits on the signal cross section and exclusion 
contours in the interpretation of the GGM SUSY breaking scenario. 



Figure 10 shows the upper limits on the GGM production cross section as well as exclusion 
contours in the squark versus gluino mass plane for the single-photon bino-like neutralino and 
the diphoton wino-like scenario. 



As further interpretation of the single-photon and diphoton results, Fig. 11 shows the 95% CL 
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Figure 10: 95% CL upper limits on the signal cross section (left) and corresponding exclusion 
contours (right) in gluino-squark mass space for the single-photon bino-like (top) and diphoton 
wino-like (bottom) scenario. The shaded uncertainty bands around the expected exclusion 
contours correspond to experimental uncertainties, while the NLO renormalization and PDF 
uncertainties of the signal cross section are indicated around the observed limit contour. 
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Figure 11: 95% CL upper limits on the signal cross section in the plane of gluino versus neu- 
tralino mass for the single-photon search in the wino-like scenario (left) and the diphoton anal- 
ysis for a bino-like neutralino (right). 



upper limits on the signal cross section for the single-photon (diphoton) wino-like (bino-like) 
scenario in the gluino-neutralino mass plane, while Fig.[l2]displays 95% CL upper limits on the 
signal cross section and exclusion contours for the single-photon (diphoton) bino-like (wino- 
like) scenarios. 



Figure 13 shows the the 95% CL upper limits on the signal cross section in the plane of wino-like 
versus bino-like gaugino mass for the diphoton analysis. 

A.2 Simplified Model Interpretation 

This section contains additional figures from the interpretation of the single and diphoton anal- 
yses in terms of Simplified Models. Acceptances for the single-photon analysis in the case of 
the W7 Simplified Model and for the diphoton analysis in the 77 SMS interpretation are shown 



in the neutralino versus gluino mass plane in Fig. 14 



The SMS results from the single-photon analysis for the 77 Simplified Model are given in 



Fig. 15 The distribution of acceptance as well as upper limits on the cross section plus ex- 
clusion contours are displayed. 

A.3 UED Interpretation 

This section contains additional figures from the interpretation of the diphoton analysis in 
terms of universal extra dimension models. The UED acceptance times efficiency is shown 
in Fig. 16 for n = 2 and 6 large ED. 



The UED cross section upper limit is shown in Fig. 17 for n = 2 large ED. The 95% CL limit for 
n = 2 is compared with the expected UED production cross sections resulting in an exclusion 
limit of 1350 GeV for an expected limit of 1340 GeV. 
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Figure 12: 95% CL upper limits on the signal cross section (left) and corresponding exclusion 
contours (right) in the plane of gluino versus neutralino mass for the single-photon bino-like 
(top) and the diphoton wino-like scenario (bottom). 
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Figure 13: 95% CL upper limits on the signal cross section in the plane of bino-like versus 
wino-like gaugino mass for the diphoton analysis. 




Figure 14: Acceptance for the single-photon analysis in the W7 Simplified Model (left) and for 
the diphoton analysis in the 77 SMS interpretation (right). 
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Figure 15: Acceptance (left) and 95% CL observed upper limit on the cross section plus exclu- 
sion contours (right) in the gluino-neutralino mass plane for the single-photon analysis inter- 
preted in the 77 Simplified Model. 
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Figure 16: UED acceptance times efficiency for n = 2 large ED (left) and for n 
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Figure 17: UED cross section upper limit for n = 2 large ED at 95% CL compared with expected 
UED production cross sections (black diagonal line). The shaded region shows the uncertainty 
due to PDFs and renormalization scale on the expected limit. 
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